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Abstract. We consider a local, renormalizable, BRST-invariant action for QCD in Coulomb gauge 
that contains auxiliary bose and fermi ghost fields. It possess a non-perturbative vacuum that spon- 
taneously breaks BRST-invariance. The vacuum condition leads to a gap equation that introduces a 
mass scale. Calculations are done to one-loop order in a perturbative expansion about this vacuum. 
They are free of the finite-T infrared divergences found by Linde and which occur in the order 
corrections to the Stefan-Boltzmann equation of state. We obtain a finite result for these corrections. 
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INTRODUCTION 

In 1980 Linde [1] showed that standard finite-temperature perturbation theory suffers 
from infrared divergences. Since then no solution has been found, although the infrared 
divergences may be avoided by introducing a magnetic mass for the gluon in an ad 
hoc manner. Some time ago, it was suggested that these divergences arise because 
the suppression of infrared gluons by the proximity of the Gribov horizon in infrared 
directions is neglected in standard perturbation theory [2]. We present an approach in 
which is based on a local action that possesses a non-perturbative vacuum. It appears 
that the infrared divergences found by Linde do not arise in this approach because the 
free propagator, in this vacuum, of 3-dimensionally transverse, would-be physical gluons 
is strongly suppressed in the infrared. 

Starting from the local action, we shall derive a gap equation that determines the 
non-perturbative vacuum. We shall solve the gap equation at high temperature to obtain 
the leading non-perturbative correction to the Stefan-Boltzmann equation of state of 
the gluon plasma. Additional results in the present approach, including propagators 
at the one-loop level, relation to the cut-off at the Gribov horizon, issues of principle 
concerning unitarity and Lorentz invariance, relation to other approaches, accord of the 
present approach with numerical studies, and references, may be found in 
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LOCAL BRST-INVARIANT ACTION 



We shall be interested in pure SU(N) gauge theory at temperature T . Finite T is described 
by a Euclidean action which, for pure SU(N) gauge theory, is of the form 

S = Jd^x^; ^ = ^YM + s^. (1) 

where ^ym = \ F^v ^^'^ F^v = d^Ay — dyA^ + gA^ x Ay. Here g is the coupling 
constant, and we use the notation for the Lie bracket (A x B)" = f'^^'^A^B'^ , where f"^^ 
are the fully anti-symmetric structure constants of the SU(N) group. The color index 
is taken in the adjoint representation, a= 1,...,A^^— 1. We shall generally suppress the 
color index, and leave summation over it implicit. Configurations are periodic in xq, 
^/i(^'''^o) = with period /3 = l/T, where T is the temperature. The 

integral over xq always extends over one cycle, / dxQ = Jq dxQ. We are in D Euclidean 
dimensions. Lower case Latin indices take values, i = 1,2,...,D— 1, while lower case 
Greek indices take values /i = 0,l,...,D— 1. 

In the BRST formulation there are, in addition to A^, a pair of Faddeev-Popov ghost 
fields c and c and a Lagrange multiplier field, b, on which the BRST operator acts 
according to 

sAn=D^c; sc = —{g/2)c X c; sc = ib; sb = 0. (2) 

It is nil-potent, = 0. Here Z)^ is the gauge-covariant derivative in the adjoint represen- 
tation, D^c = d^c + gA^ X c. 

The choice of the gauge-fixing density in ([T]) is the choice of gauge. Physics is 
independent of ^, provided that it results in a well-defined calculational scheme. For 
finite T, this is a crucial proviso, because the standard gauge choice leads to infrared 
divergences fT]. The standard Coulomb gauge is defined by ^ = ^coui = dtcAi, with 
= idibAi — djcDiC. The Lagrange-multiplier field b imposes the Coulomb gauge 
condition djAi = 0. To avoid the infrared divergences of the standard gauge choice we 
introduce an additional i'-exact term s^aux that involves a quartet of auxiliary ghost fields 
on which s acts trivially 

5<^);* = <; 5< = 0; sa)f = r^'; = 0. (3) 

The fields and are a pair of bose ghosts, while co^^ and O)^ are fermi ghost 
and anti-ghost. The indices a and b label components in the adjoint representation of 
the global gauge group, a,Z7 = 1, ...,A^^ — 1, and /i is a Lorentz index. We take for the 
gauge-fixing term 

^ = Ud + = diC^Af + d,cdf{Di(^^r\ (4) 

where we stipulate that the gauge covariant derivative acts on the first color index only, 
[Di^^^r' = ,9,0^* +^/-'^A^-(^)f etc. 



NON-PERTURBATIVE VACUUM 



Maggiore-Schaden shift 

We make a change of variable whereby the bose ghosts are translated by a c-number 
term linear in the spatial coordinate 

< W = < W - -^m; W = W + r'^' ^m- (5) 

Note that and cp designate different fields. Here 7 is a parameter with dimensions of 
(mass)^ that will be determined by the condition ^ = 0, where W is the free energy. 
We also translate b and c by compensating terms, 

which are chosen to cancel explicit .Jc-dependence in the new action. The BRST operator 
s acts on the new variables according to 

^'P? = «?; 5d)«* = 9;;* + //2 sc"^ = ib'''^; sb""^ = 0. (7) 

Despite the ^-dependent shift, remarkably, neither the gauge-density dU nor the La- 
grangian density acquire any explicit jc-dependence when expressed in terms of the new 
variables, 

^ = d.t^'A^ + (A(Pm)"' - y''^ (D^cb^r. (8) 

As before, it is understood that the gauge-covariant derivative acts on the first index 
only, {Di&^Y^ = diCJ^^ + gf'^^A'ioi}^'' . After the shift, the Lagrangian density is given 
by ^ = + st, , where 

= idib'^^A'^ - <9,c-*^XAc)" + ^,-<(A(Pm)"' - {D.co^r' + {8D,c x (p^)"' ] 

+7I/2 (D.cpiY'" - 7I/2 [ {D,q>ir" + {gDiC X w,-)"" ] - 7 - 1)(^ - 1), (9) 

and (D,c x 9^)"* = f"-^'^[Dicy ^'^^^ acts on the first color index, etc. For purposes of 
expansion in powers of g we shall change indepdendent parameter from yiom according 
to 7^/2 = (^2N)U2g ' where m has dimensions of mass, and m is taken to be of order g^. 

Gap equation 

Henceforth we shall be concerned with the action 5, regarded as a function of 
the new fields. The partition function is given by Z = / J4>exp(— 5), where ^» = 
(A^, c, c, b, (p, ^, (0, Cb) is the set of all (new) fields, and J4> represents integration over 
them. For simplicity we have written c and b instead of c* and b*. The field cp^^ is real 
while (p^'' is pure imaginary. The classical vacuum occurs where all these fields vanish 
4> = (A^,c,c, Z?, (p, 9, £0, cb) = 0. Finally the value of 7 is determined by the condition 



that the free energy W = InZ be stationary, = 0, or i^-^^j = 0. There is a non- 

perturbative vacuum if this equation has a solution with 7 7^ 0. We do not require that W 
be a maximum because there are unphysical fields present. The last equation reads 

- ^>ir - isD,c X cd,r) = 2//2(A^2 _ 1) (z) _ 1). (10) 

The second term vanishes, {{gDiC x cbi)"") = 0, because there is no cO) term in the action. 
Moreover the new Lagrangian density ^ is invariant under space-time translation of the 
new fields, ^{x) 4>(x + a), and the vacuum just found, at 4> = 0, is also. Translation 
invariance implies that the terms di(p and di^ do not contribute to (flOl) . and we obtain 

{r'^-AUcp-cpyn = ^ {D-i){N^-i). (11) 

This gap equation determines m = m{g,T). Invariance under scale transformation is 
spontaneously broken for m 7^ 0. 



Spontaneous breaking of BRST symmetry 

The quantity is the analog of the vacuum expectation-value v of the Higgs field 
4> that appears in spontaneous symmetry breaking of global gauge invariance. For in 
the Higgs mechanism one makes the translation = 4>*'^ -|- vd^, and the vacuum 
expectation-value v is determined by the condition that the free-energy be stationary 
with respect to v, ^ = 0. Similarly, trp- is determined by the condition ^ = 0. In 
the present case, BRST invariance is spontaneously broken rather than global gauge 
invariance, because the expectation-value of i'-exact quantities is non-zero, for example, 

( s(bf ) = {q>f + ) = x^ ^ 0. (12) 

As in the Higgs case, the spontaneously broken theory inherits renormalizability 
from the unbroken theory. But (fT2l) shows that we cannot identify observables with 
equivalence classes of j'-invariant objects, modulo i'-exact quantities, as in the standard 
BRST approach. However, as shown in [3], the present method is formally equivalent to 
the canonical formulation of Coulomb gauge, with a cut-off at the Gribov horizon. This 
allows us to identify observables, such as the energy-momentum tensor T^y^ with the 
corresponding quantities in the canonical formulation so, for example, T^y = F^X^v ~ 



FREE PROPAGATORS 

We now develop a perturbative expansion about the new, non-perturbative, vacuum. 
For this purpose we treat m as an independent parameter of order g*^, and calculate 
perturbatively all one-particle irreducible graphs, including the gap equation, to a given 



order in g. Then m = m{g, T) is determined by solving the the gap equation (fTTI) non- 
perturbatively. 

The first step is to expand the action in powers of g, S = 5_2 + 5o + ... . The leading 
term is of order g^^, 

S-2^-:^{N^-l){D-l)L'P, (13) 

where the spatial quantization volume isV = L?, and the time extent /3 = T ^ Although 
5 2 is independent of the fields, and does not contribute to the propagators, it should 
not be ignored because, when the gap equation is solved for m = m{T), it gives a T- 
dependent contribution to the free energy. The terms in the action of order are all 
quadratic in the fields 

So= Jd^x \^d^A^y-dyA'l,)^ + idib''A^-d,c''dic'' 



,2 



(14) 



and determine the "free" propagators. The term with coefficient rrp- causes a mixing of 
the zero-order transverse gluon and bose-ghost propagators. 

To calculate the free propagators, we define the field, = ^2Ny/2 f"^'^i^j" ~ 9)°)' 

that mixes with . The orthogonal component j^^j^yjif^^^'i^j" + 9)") other com- 
ponents of (p and 9 do not mix with A/. The free propagators are given by 

DA^AjiKko) = DA,y,j{Kko) = D^,^.{k,ko) = '\ (15) 

where A = (k^ -|- kQ)k^ + m^, and Pij = 5ij — kikj is the transverse projector. Here ko — 
Inn/P are the Matsubara frequencies, where n is any integer, and we have suppressed 
the trivial color factor 5'^^'. In terms of the variable xjfj, the gap equation (fTTI) reads 

-/ {A%0)rjm = ^ {D-l){N^-l). (16) 

GAP EQUATION IN ONE-LOOP APPROXIMATION 

When the left-hand side of the gap equation is evaluated to zeroth order in g, using the 
mixed propagator (fT5l) , it reads 

d^-^k y D-2 

(2;r)0-i |;(k2 + yt2)k2 + m4 " A^' ^ ^ 

where we used /'//(k) = D — 2. The sum over Matsubara frequencies yields 

d^-^k (D-2) / 2 ^_D-l 

(2;r)^-i 2k2£ V ^ expipE) - 1) Ng^ ' ^ ^ 



which holds in D Euclidean space-time dimensions. The first term in parentheses gives 
the result at T = 0, and the second term is a Planck-type finite-temperature correction. 
For D < A the integral is convergent. We take the limit Z) — > 4. The first term in 
parentheses has the limiting form 



d^-^k (Z)-2) 1 



(2;r)o-i 2k2£ An^ 



1 , 



-+lnM-^ 



(19) 



where £ = (4 — D)/2, and /i = is, in general, a temperature-dependent renormal- 
ization mass. There is a pole at D = 4. We subtract the pole term, and the gap equation 
reads 

1 rdx 1 _ 'iTt^ 

I'^KmJ^jQ 'V exp{ml5u)-l~ Ng^i^y ^ 

where u = {x^ + ^) ^Z^. 

We now specialize to high temperature T, and take the renormalization mass to be 
II = T. This simplifies the high-T case because of asymptotic freedom. The running 

coupling g{T) is small, and is given approximately by -^jj^^ = ^^-^ In (j^^^ where 

Aj^ is the MS physical mass scale. At high temperature, the first term in (|20l) may be 
neglected, and the second term has the limit 

dx 1 1 /-""^^ _ TtT ^21) 

u exp(/3mM) — 1 /3m Jo X^ + l 2^/2 m' 

so the gap equation at high T simplifies to = tJ'^lj) ' ^^^^ solution, 

'"^^'^) = ^37^ s\t) t r ^ oo. (22) 

Thus, in the high-temperature limit, m{T) is proportional to the magnetic mass g^{T)T . 

FREE ENERGY 

To order g^, the free energy W = InZ is given by 

expW = J d^ exp{-S-2-So), (23) 

where given in (fT3] ). is field-independent and of order g~'^, while Sq, given in (fT4l) . 
is quadratic in the fields. We obtain W = W-2 + Wq where, for Z) — 1 = 3, W-2 = —S-2 = 

3m'' fj^T2 ^\T3 



irrp2 {N — l)L^/3, and expWb = / J4>exp(— 5o). The evaluation of the free energy Wq 
for the quadratic action 5o is straightforward, keeping in mind the A-cp and A-(p mixing, 
with the result 

(A^2_i)y/3 p (fc4_^4) 
^° ^ 3;r2 io £[exp(£/3)-l] ' ^'^^ 



where E = E{k) = {k^ + ^Y^^- We add the term W-2 and obtain to order the free 
energy per unit volume, w = W /V, 

EQUATION OF STATE AT HIGH TEMPERATURE 

We now evaluate w in the high-temperature limit, where g{T) is small, and our expansion 
should be reliable. We insert the solution, m = ~j§jTY^' S^P ^^l^ation and obtain 

where K(n) s V u s + 4) "\ and „ ^ (f )" = (|^)' is a small 

parameter. With neglect of higher order terms, one obtains K{r}) = ^ — ^^/i , which 
gives 



w ■ 



This is the leading correction to the Stefan-Boltzmann limit from the non-perturbative 
vacuum. The equation of state of the gluon plasma follows from the thermodynamic 
formulas for the energy per unit volume and pressure, ^ = P = and entropy 

per unit volume, s = To calculate the energy density, we use — = = /3- 

function = 0{g^), which is of higher order. We thus obtain for the energy density and 
pressure at high temperature, e = 3p = w{T)T, and s = |w(r). 

We have obtained the leading contribution that comes from the non-perturbative 
vacuum. Numerically it is a small correction, whereas the correction of order g^ is 
divergent when calculated with the perturbative vacuum To this must be added the 
perturbative contributions, including resumations, that have been calculated at m = 0, 
and that are of lower order in g [5]. Since standard, resummed perturbation theory 
diverges at order g^, which is precisely the order of the correction we have found, the 
result obtained here is consistent with standard perturbative calculations. 
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